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I discuss the size distribution N{S) of avalanches occurring at the yielding transition of mean 
held (i.e., Hebraud-Lequeux) models of amorphous solids. The size distribution follows a power 
law dependence of the form: N{S) ~ S~'^. However (contrary to what is found in its depinning 
counterpart) the value of r depends on details of the dynamic protocol used. For random triggering 
of avalanches I recover the r = 3/2 exponent typical of mean held models, which in particular is 
valid for the depinning case. However, for the physically relevant case of external loading through 
a quasistatic increase of applied strain, a smaller exponent (close to 1) is obtained. This result 
is rationalized by mapping the problem to an ehective random walk in the presence of a moving 
absorbing boundary. 

PACS numbers: 


I. INTRODUCTION 

The yielding transition (at zero temperature) of an 
amorphous solid material occurs when the externally ap¬ 
plied shear stress a overpasses a critical value CTc- For 
a < Gc the system is blocked and remains in this state 
indefinitely. As cr > CTc the system enters a plastic flow 
regime in which the strain in the system increases linearly 
with time. This is referred to as a plastic flow regime. 
The transition at tr = CTc is the yielding transition. At Cc 
the system is in a critical state, and the dynamics pro¬ 
ceeds via avalanches of all sizes, characterized by its size 
distribution M{S). 

The previous scenario has been confirmed by a vari¬ 
ety of experimental measurements A” well as analyti¬ 
cal and numerical techniquesFrom all these sources, 
the idea has emerged that plasticity proceeds through 
the destabilization and rearrangement of discrete units 
at given spatial positions (sometimes referred to as shear 
transformation zones) d^i^° Rearrangement of a unit at a 
given spatial position has an effect (through direct elas¬ 
tic interaction) on the stability of other units at different 
locations. It is thus clear that the form of the elastic in¬ 
teractions between different regions of the material play a 
central role in the problem. In this respect, it is crucial to 
realize that the elastic interaction kernel has alternating 
signs in different spatial regions^^J The destabilization of 
a given site can thus have both a destabilizing or stabi¬ 
lizing effect on some other site. 

The previous description shares many features with the 
well known case of depinning transitions of elastic man¬ 
ifolds on disordered energy landscapesHowever, a 
crucial difference steps from the fact that for depinning, 
destabilization of a given site always produces a desta¬ 
bilizing effect on any other site, namely the interaction 
kernel is positive definite. The mean field description of 
the depinning transition is well studied, and is known to 
produce a critical size distribution of avalanches of the 
form M{S) ^ It has to be mentioned that this 

value of r stands for (at least) two quite different proto¬ 


cols used to trigger avalanches in the system. In one case 
[random loading) one site is chosen at random and the 
local stress is increased until the site destabilizes. In the 
second case [uniform loading) the stress on the system 
is increased uniformly over all sites, until the first site 
becomes destabilized. In both cases r = 3/2 is found for 
depinning. 

I will show that in the case of mean field plastic yield¬ 
ing models, the value of r that is obtained is dependent 
on the loading mechanism applied. While for random 
destabilization the value r = 3/2 is obtained, uniform 
loading produces a smaller value of r, close to 1. The 
reason of this difference is the much larger stress that is 
added to the system in the uniform loading protocol com¬ 
pared to the random loading case. I present the model to 
be studied in the next Section. An analytical treatment 
of the problem is given in Section III, and results of nu¬ 
merical simulations are presented in Section IV. Section 
V contains the conclusions. 


II. THE MODEL 

Hebraud and Lequeux^i considered a mean field model 
of the plastic yielding transition. The model considers a 
real valued variable yi, representing the stress at each 
spatial position of the system.?- The coupling among yi 
variables is mean field like, i.e., any change in yi affects 
equally all other y’s. I discuss here a small variation of 
the original HL model, in the limit of quasistatic loading. 

I consider N stresses yi, i = which define the 

instantaneous configuration of the system. I will refer to 
the N variables i alternatively as ” sites”, or ” particles”. 
The configuration is stable as long as all yi are lower 
than a plastic threshold value, taken for simplicity to be 
uniform, of value 1. It will be more convenient to describe 
the model in terms of the variables Xi = 1 — yi. Each Xi 
represents the amount to additional stress necessary to 
destabilize site i. The dynamics proceeds according to 
the following rules. 
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(a) On a stable configuration, an instability is pro¬ 
duced (either by random or uniform loading), and an 
avalanche starts. Suppose i is the site that becomes un¬ 
stable. 

(b) Xi is increased a quantity (1 -I- k)u, where m is a 
random quantity, of order 1 (below I take u to be ex¬ 
ponentially distributed: p(u) = exp(—u/u), u > 0 ), 
and fc is a non-conservation parameter (in such a way 
that criticality is expected as fc —>• 0 ). 

(c) All Xj are decreased a quantity ■q/'/N—u/N, where 
77 is a Gaussian variable, with zero mean and fixed stan¬ 
dard deviation 

(d) : Steps (b-c) are repeated over all destabilized 
sites^ until all sites become stable (this defines the end 
of the avalanche). Then I return to (a). 

The instantaneous stress a in the model is defined as 
cr = — Xi). The size S of an avalanche is con¬ 

veniently defined as the sum of all u values generated 
during the avalanche. Note that according to the previ¬ 
ous rules, an avalanche of size S generates a reduction 
of stress in the system of value kS/N. As fc —>■ 0 the 
average size of avalanches goes to infinity as l/k. This 
corresponds to a broad distribution of avalanches sizes in 
this limit, which corresponds to criticality. For finite k, 
there is typically a maximum size Smax of the avalanches 
that are generated. 

The previous rules indicate that each site reaching the 
a: = 0 border is reinserted at x = {1 + k)u, then a: = 0 
is an absorbing border. Each unstable site also produces 
random kicks (of intensity q/'/N) plus a systematic de¬ 
crease (of value —u/N) on all other sites. In terms of a 
fictitious time that counts the number of sites that reach 
the a; = 0 border, these two terms can be described as 
diffusive, and convective, respectively. 

In the thermodynamics (large N) limit, the equilib¬ 
rium distribution of x values is characterized by a func¬ 
tion P(a;), such that NP{x)Ax gives the number of sites 
with X within the small interval Ax. The form of P{x) 
is determined as a balance between the sites that reach 
a: = 0 and are thus redistributed, and the change in x that 
each site reaching a; = 0 produces. In an equilibrium sit¬ 
uation the total change AP produced by a single particle 
being destabilized and reinserted must vanish, and this 
means that P{x) satisfies the equation (for A: —>■ 0) 

( 1 ) 

with D = ql/2 of step (c) above, and P{x) = 0 for a: < 0. 
The three terms of the r.h.s. represent respectively the 
effect of diffusive and convective variations of x (the two 
terms in the change of Xj at step (c)), and the effect of 
reinserted sites (step (b)). The amplitude of the rein¬ 
serted term, namely DP'{0) + P{0) takes into account 
that sites can reach the a; = 0 border due to the diffu¬ 
sive, or convective evolution. 

There is an important difference in the equilibrium 
form of P{x) when D = 0, or D > 0. In the first case 


(in which the previous model describes depinning) P{x) 
is immediately obtained from o as 

rOO 

P{x) = -P(O) / p{x)dx (2) 

J X 

In particular, P(x) goes to a non-zero value when a; —>■ 0. 
For £) > 0, due to the diffusive term and the absorb¬ 
ing boundary condition, P{x) goes to zero at a; = 0. 
Moreover, the particle rate across the a; = 0 border due 
to the diffusive term is proportional to dP/dx\x=Q, and 
for this quantity to be finite, the form of P{x) must be 
P{x) ^ Ax for small x. This is a mean field results that 
does not hold for short range interactions, where a form 
P(x) ^ x^ with a non-trivial 0 is expected. 

For fc —>■ 0, P(x) is independent of the form of the 
destabilizing mechanism (random or uniform loading) 
since at criticality the stress redistribution caused by di¬ 
rect external loading is negligible compared with reac¬ 
commodations during avalanches (this is not true for fi¬ 
nite k, see the Appendix). The form of P{x) for different 
values of D are presented in Section IV. 


III. RESULTS FOR AVALANCHE STATISTICS 

In order to calculate the avalanche statistics, the con¬ 
tinuous description provided by the P{x) function is not 
sufficient. In fact, what this distribution tells is that each 
single site reaching a; = 0, produces changes in all other 
X such that on average an additional site reaches a: = 0, 
i.e., a conservative critical situation is reached. However, 
in order to describe avalanche distribution we must con¬ 
sider fluctuations. 


A. Random loading 

It is convenient at this point to describe the dynamics 
of the system in term of a discrete ’’time unit” corre¬ 
sponding to one step of points (b) and (c) of the dynam¬ 
ical protocol described above. Namely, at each discrete 
time, one kick is given to every particle according to (c), 
and one destabilized particle is reinserted according to 
(b). Let us call Ui the cumulative number of destabilized 
particles until time t = i. The condition for an avalanche 
started at t = 0 to survive until time t = i is that Ui > i. 
The avalanche stops at the moment in which Ui < i. 

Since the sites carrying different values of x, once 
destabilized, are reinserted through a random process 
(implied by the random values of the u parameter) the 
exact values of xt are locally uncorrelated. Due to this 
uncorrelated behavior, the stochastic process ni corre¬ 
sponds to a Poisson process with a unitary rate per time 
step. Dehning the stochastic process as mi = ni — i, 
we see that the avalanche lasts until the hrst time iq in 
which mig = 0. Note that as one single site is reinserted 
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at each time step, io measures also the avalanche size, 
i.e., S = lo¬ 
in the large i limit rrii is a symmetric random walk, 
with probability distribution p{mi) given by 


= Tin 

and the previous description corresponds to its survival 
probability in the presence of an absorbing boundary at 
the origin. This surviving probability goes as The 

size distribution of avalanches is given by the probability 
density of the random walk being absorbed at time t, 
which gives M{S) ^ , i.e., r = 3/2. I emphasize 

that the previous method shows this result is valid for 
random loading, both in the case tjq = 0 (depinning), and 
r]o > 0 (plastic yielding) as it only depends on the fact 
that the number of sites reaching the instability threshold 
is an uncorrelated stochastic process with constant rate. 


B. Uniform loading 

In the uniform loading case, in order to start an 
avalanche we look for the smallest x value in the system, 
say Xmin, and add this additional load to the system. 
This is equivalent to say that all cc’s in the system are re¬ 
duced in the quantity Xmin- The question is if this shift¬ 
ing of the P{x) distribution may bring some observable 
effect. We will see that the answer depends on whether 
we are considering the case 770 = 0 or t^q > 0 . 

For r]o = 0 (depinning case) the form of P{x) tends 
to a constant as a: —>■ 0, i.e, P{x) ^ P{0) for small 
X. When shifting the distribution by Xmin to start the 
avalanche, we modify P{x) to a new P{x) given by 
P{x) = P{x -\- Xmin) P(0) (to leading order). This 
means that the form of P(x) is not greatly affected by 
the Xmin shift. In fact, results of numerical simulations 
described below show that both loading protocols pro¬ 
duce the same avalanche distribution A/’(5') ^ in 

the r]o =0 case. 

However, the situation is different in the case rjQ > 0. 
Now P{x) ^ Ax for small x, and a shift of the distri¬ 
bution by Xmin transforms it to P{x) = P{x + Xmin) ^ 
Ax + Axmin, namely the probability distribution gets a 
constant correction near a; = 0 which qualitatively mod¬ 
ifies its form. This produces an important change in the 
avalanche size distribution as I will show now. 

Suppose we have the distribution P{x) = Ax + Aq. 
The linear part generates a flux of particles through the 
X = Q border characterized by the Poissonian stochastic 
process n^. The average rate of this process is easily 
seen to be equal to AD, and since we have one particle 
arriving per unit time, we conclude that A = 1/D. The 
constant term Aq produces the arrival at a; = 0 of some 
additional particles, that contribute to the process. Let 
di be the number of additional particles arrived at a; = 0 


at or before time i, originated in the additional constant 
term Aq. Now the stochastic process to be considered 
in order to determine the size of the avalanches is fhi = 
rii — i + di, i.e, the avalanche survives as long as fhi > 
0. We need to characterize the contribution di from Aq. 
Due again to the fact that sites that contribute to the 
constant Aq are totally uncorrelated, di is a (cumulative) 
Poissonian process, however now its average rate depends 
on time. The time dependence of the average rate can 
be calculated simply noticing that it gives the average 
number of particles absorbed at a; = 0 , starting from a 
distribution with constant value Aq for all a: > 0. The 
average number di can be obtained by direct integration 
of the diffusion equation, and the result is 



The condition for an avalanche to survive until time i 
is then fhi = rrii ~\- 2 Aq^JD i/n > 0. In other words, 
the original condition Wi > 0 is now replaced by rrii > 
- 2 H 0 /DI 7 TT, i.e., the absorbing boundary that was orig¬ 
inally at i = 0 now retracts in time, as ^ 

The effect of such a moving absorbing boundary con¬ 
dition on a random walk has been studied in the litera¬ 
ture (see for instance Refs^“— ). The important point 
that makes the problem analytically solvable is that the 
position of the absorbing boundary has the same depen¬ 
dence in time that the internal diffusive dynamics, and 
this allows to transform the problem to a single ordinary 
differential equation.— The survival probability is thus 
observed to remain a power law in time, but with a non- 
universal exponent. The value of r that depends on the 
coefficient R = 2Aq ^D/-k of the square root time behav¬ 
ior in Eq. |4] (assuming a normalized underlying random 
walk, as in Eq. [3]). It is found that r = 3/2 for i? = 0, 
and T —>■ 1 when i? —>■ 00 . We see the full form of this 
dependence in Fig. [TJ 

In order to calculate the actual value of R, we must 
determine what value of Aq is appropriate in our case. 
It turns out that we do not have a unique value of Aq, 
but a continuous distribution. In fact, given a starting 
distribution P{x) = Ax, the lowest x value, which defines 
Xmin is itself a random variable, with distribution 


p{Xmin) = Axmin exp ^ (5) 

We can write Xmin = Xq/ VA, where xq follows the nor¬ 
malized distribution 

p(xo) = a;oexp (6) 

The values of Aq are obtained from Aq = Axmin = 
VAxq, giving for R the value (remember that A = l/D): 



( 7 ) 
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FIG. 1: The r{R) function obtained from the mapping to the 
survival probability of a random walk of variance = t, in 
the presence of an absorbing wall that recedes as -Ry/i. I also 
show the distribution of R values according to its probability 
distribution p{R), and the value of t ~ 1.09 obtained for its 
average value. 


We can estimate a typical value of r as the one corre¬ 
sponding to xo equal to its average value (xg = y/iFpl). 
This gives R = y/2, that entered in Fig. [Hgives r ~ 1.09. 
However we must notice that different avalanches are pro¬ 
duced with different xg values, and the final distribution 
of avalanche sizes has to be calculated as a composition 
from avalanches taken from distributions with different 
T values according to 


Mtot{S) = / dTp(T)N {S{t)) ( 8 ) 

JSmin 

where p ( t ) = p{xg)dxg/dr (taken from data in Fig. [T|) is 
the probability of different values of r, and N{S{t)) = 
CtS~'^ , with Cr a normalizing constant. Taking the 
minimum size of avalanches to be Smin = 1, it results 
Cr = (r — 1), and a numerical evaluation of the expres¬ 
sion (|8l) produces a non-perfect power law, that can be 
characterized by an effective, local exponent defined 
as 


'^eff — 


dlog[iVTOT(«S')] 

dlog^ 


(9) 


The result of this evaluation can be seen in Fig. [2j As 
we move to larger values of 5, T^a is more dominated 
by the terms with the smallest r in Eq. (|S]), and 
decreases. It has to be noticed however, that for this 
effect to be appreciable we must really move to extremely 
large values of S. As a rule of thumb, we can say that 
we expect a t between 1.1 and 1.2 to be observed. 



FIG. 2: The local exponent Tetf as a function of the avalanche 
size S. The value of Tett is a very weakly decreasing function 
of S. 


IV. NUMERICAL SIMULATIONS 

I will present here results of numerical simulation of 
the model described in Section II, that will confirm the 
scenario described in the previous Section. In all the 
simulations presented the values of u are taken from an 
exponential distribution with mean value of 1, namely 
p{u) = exp(—u). In this case, in the thermodynamic 
limit the equilibrium form of P{x) can be readily found 
by integrating Eq. (TJ and is given (for fc —>■ 0) by 

The form of P{x) for different values of rjg at criticality 
{k —^ 0) is shown in Fig. [S] Note that it must be < 1 
for the solution to exist. As it was explained before, 
the main feature of these curves is their behavior near 
X = 0. We see that that P{x) ^ x° for r]g = 0, whereas 
P{x) ^ X^ for T]g ^ 0. 

Next, I present results for the avalanche size distribu¬ 
tions. Results in Fig. H] correspond to the depinning 
model {rjg = 0). They were obtained in systems with 
N = 10®. It is seen that the size distribution AfiS) 
is a power law Af{S) ~ 5“®/^ that is cut-off at large 
avalanche size at some characteristic typical maximum 
size Smax, that increases as k decreases. We can obtain 
the increase law of Smax with k as follows. An avalanche 
of size S generates a reduction in the systems stress of 
value Sk. This quantity is (on average) compensated by 
the stress increase in the loading stage. For the r]g = 0 
case, this increase is given on average by u, which is equal 
to one in our case. So we get S = k~^. Also, S is related 
to Smax through S 5'^“^, which gives Smax ^ k~~, 
and in the present case with r = 3/2, we get Smax ^ k~^, 
that is well satisfied by the results in Fig. 0] Note that 
the results for rjg = 0 show no difference for the two differ¬ 
ent driving protocols, namely random or uniform loading, 
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FIG. 3: Form of the function P{x) for different values of t]q, in 
the fc —>■ 0 limit. The left panel contains results of simulations 
with N = 10000. The right panel shows the analytical limit 
for N ^ oo. The crucial characteristic of these curves is the 
behavior near x = 0. Note that P{x) ~ x° for rjo = 0, whereas 
P{x) ~ x^ for Tjo / 0. 



FIG. 5: Same as previous figure but for » 7 o = 1 in the random 
loading case. The results are seen to be similar to the case of 
depinning (770 = 0), with an avalanche size distribution with 
an exponent r = 3/2, and a large size cut off Smax ~ 


I notice the rather strong overshoot of the size distribu¬ 
tions near Smax, that I analyze in the Appendix. 


according to the arguments given in the previous Section. 



FIG. 4: Simulations in a model with A = 10®, and 770 = 0 (i.e., 
depinning case). Avalanche size distribution for decreasing 
values of k. Results for the two different loading protocols 
(uniform and random) are shown, and they are seen to be 
equivalent. The distribution observed has the form M{S) ~ 
5“®^^ exp{—S/Smax). The dependence of Smax on k is shown 
in the inset, and is Smax ~ k~^. 

Now we move to the analysis of the plastic yielding 
case (? 7 o 7 ^ 0). In this case, the two loading protocols 
generate different results. For random loading we find a 
power law distribution with r = 3/2 (Fig. [5]), equivalent 
to that found in the depinning case. 

In the case of uniform loading the results are different, 
as expected from the analysis of the previous section. In 
Fig. ini we see the decay of number of avalanches with 
size. On the range of sizes displayed in Fig. [51 a decay¬ 
ing exponent r ~ 1.1 can be estimated as k is decreased. 



FIG. 6 : Simulations for 770 = 1, A = 10® in the uniform 
loading case. It is seen that as k decreases, a power law with 
a value smaller than the standard r = 3/2 value is observed. 
It is also apparent the existence of a rather strong overshoot 
at the range on the largest avalanches observed. 


To gain some additional insight on the results in Fig. 
[ 6 l I present some additional analysis concerning the 
avalanche size and the values of Xmin at which they are 
triggered. In Fig. [7] we see the distribution of values of 
Xmin obtained along a simulation with N = 10®, k = 0.1. 
We see that the distribution follows very closely the dis¬ 
tribution given by Eq. ®, which in turns originates in 
the linear distribution of x values near x = 0. 

In addition, if avalanches are independent, and they 
follow the results from the theory of the previous sec¬ 
tion, we should expect a different decaying exponent for 
avalanches if they are classified according to the value of 
Xmin at which they were triggered. In order to check for 
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FIG. 7: (a) The distribution of values of Xmin in a simulation 
with N = 10® and rjo ~ k = 0.1. The continuous line is the 
analytical expectation given by Eq. 

this, in the inset to Fig. |8] I show partial size distribu¬ 
tions of avalanches restricted to given intervals of values 
of Xmin that trigger them. If we concentrate in the low 
size part of this plots, the results follow a trend compati¬ 
ble with the theoretical analysis: avalanches started with 
smaller Xmin have a distribution with a decaying expo¬ 
nent close to 3/2, while those with larger Xmin produce 
a distribution with a slower decaying exponent. There 
seems to be however systematically lower values in the 
simulations compared with the theoretical results. 



’min 


FIG. 8: The r exponent of partial distribution of avalanches 
in the simulation of the previous figure, classified according 
to the value of Xmin at which avalanches are triggered. In the 
inset we see an example of the distributions and the low size 
ranges (thick straight lines) used to extract the r exponent. 
The three curves in the inset correspond to the three points 
in the main plot indicated by arrows. 

There is an additional effect in the uniform loading 
plastic yielding case that deserves to be mentioned. For 
depinning, or random loading plastic yielding, I have 
argued that the mean size of the avalanches must sat¬ 
isfy the relation Sk since the average external load 


between avalanches is order 1. In the case of uniform 
loading in plastic yielding, the average loading between 
avalanches is order Nxmin ~ y/V. So in this case we 
have Sk ^ \//V, or S ^ ^/N/k. Introducing the size 
of the largest avalanches in the system Smax we get 
Smax ^ {y/N This indicates that the size of 
the largest avalanches in the system are not completely 
determined by the value of k. There is also a direct de¬ 
pendence on N. In particular, at any fixed value of /c, 
Smax becomes arbitrarily large as IV —>■ oo. 

V. CONCLUSIONS 

In this paper I have studied the avalanche size dis¬ 
tributions in mean field models of the plastic depinning 
transition. The main conclusion is that contrary to the 
case of elastic depinning where the mean field exponent 
is well defined and given by r = 3/2, here the result de¬ 
pends on details of the driving mechanism. I have shown 
that for plastic depinning, the r = 3/2 exponent is recov¬ 
ered in the case of random loading, where each avalanche 
is started by making unstable a random site. However, 
the more standard loading mechanism corresponding to 
increase uniformly the stress in all the system until an 
instability is obtained, produces a significantly lower ex¬ 
ponent T ~ 1.1 — 1.2. The origin of this effect was elu¬ 
cidated by mapping the present problem to the survival 
probability of a random walk in the presence of a moving 
absorbing boundary. In addition, this mapping provides 
insight into the size crossover for cases in which the sys¬ 
tem is not right at criticality. 

The present results call the attention on the precise 
destabilization mechanism used in different models of 
plastic depinning. It is likely that the effects discussed 
here remain in more realistic cases (not mean field) where 
a more realistic symmetry and interaction range of the 
elastic kernel is considered. The kind of variation of the 
T exponent in those cases remains to be elucidated. 
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Appendix A: The approach to criticality, and the 
large size overshoot 

It has already been noticed (in connection with data 
in Fig. El) the important size effect observed at large 
avalanche sizes, where an excess number of avalanches 
(compared with the results for random loading. Fig. E]) 
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is observed. In order to study this effect in more detail, 
it is not enough to consider the critical {k = 0) case, but 
it is necessary to analyze the form in which the system 
evolves when a finite k is reduced towards 0. 

I restrict here to the plastic depinning case (770 0). 

The modifications in the balance equation o produced 
by a finite k are different in the cases of uniform, or ran¬ 
dom loading. In the case of random loading, a finite 
k produces (according to the rules of Section II) that 
the reinsertion does not occur at the position u, but at 
u{l+k). This alters the argument oip and the normaliza¬ 
tion of the last term in ©■ In addition, for finite k, the 
average size of the avalanche =llk is finite, and we must 
take into account that 1 every l/k particles is destabilized 
by the random mechanism that initiates the avalanche. 
As this random selection is made on the actual distri¬ 
bution P{x), this produces a negative term —kP/N per 
destabilized particle in the balance equation, that finally 
reads 


0 = 


D d^P 
N dx'^ 


1 dP 
N dx 


kP DP'{0) 


N A^(l 


k) 


p(x/(l -fi k)) 


(Al) 

In the case of uniform loading, we still have the change 
M —>■ m( 1 -f fc) in the reinsertion term, but now the uni¬ 
form loading mechanism produces a shift of P in a quan¬ 
tity Xmin every S destabilized particles. This gives a 
correction SP in the balance equation that is equal to 
6P = {dP/dx)dx = {dP/dx){xmin/S) = {dP/dx){k/N) 
where the last equality follows (on average) from the bal¬ 
ance between the stress increase during loading (xmin), 
and the stress decrease during avalanche (kS/N). Fi¬ 
nally, the balance equation for uniform loading takes the 
form 


D d^P l + kdP 
N dx"^ N dx 


DP'{0) 
N{l + k) 


p{x/{l + k)) 


(A2) 


For our choice p{x) = exp(—x) the solutions to both 
equations are 


P{x) 


_ x(l+Vl+ikD) 

e 1 +'“ — e 


-{ -L k _ 21? 

^ l + y/l+AkD 


for random loading, and 


(A3) 


P{x) 


_ X a: (1 + fc) 

-T+v _ e 13— 

1 + ^-m 


(A4) 


for uniform loading. 

In the limit N ^ 00 , the distribution of avalanches is 
controlled by the form of these expressions near a; = 0. 
Close to criticality, we must search for the leading terms 
in an expansion in powers of k. For random loading, lin¬ 
earizing Eq. (IA3I) near a; = 0, and writing the coefficient 
to first order in k we obtain 


P(x) = ^(l-fc(l-P)) (A5) 

where we see that as fc —>■ 0, P{x) tends to the critical 
form P{x) = x/D. The finite value of k provides a cut 
off in the size distribution of avalanches. In fact, in the 
mapping to the random walk problem, the linear in k 
term of the previous expression corresponds to an ab¬ 
sorbing moving wall at rrii = ki/D (in the notation of 
Section III). This problem can be analytically solved in 
the large avalanche size limit, and provides an avalanche 
distribution given by 

N{,S) - exp(-S'fcV^^) (A6) 

This gives a size cutoff Smax that increases a I/fc^ as 
fc — 7> 0, and results that fit nicely the curves in Fig. [5l 

For uniform loading instead, it is easily verified that 
Eq. (IA4I) is independent of k to first order in x^ This 
means that in order to obtain an avalanche cutoff we 
must expand to second order in x. The result (to first 
order in k in the coefficient of x'^) is 

Pix) = ^-^{D + l + k{l-D))x^ (A7) 

This form of P{x) allows to explain in detail the be¬ 
havior observed in Fig. |6]for the uniform loading case. In 
fact, when k = 0 we are at criticality, and we can describe 
the avalanche distribution through the model of Section 
III, of the survival of a random walk, where the loading to 
Xmin corresponds to the presence of a retracting absorb¬ 
ing wall at ^ —XminVt- The linear in k term in Eq. (IA7I) 
corresponds to an additional current of particles through 
the absorbing wall that has to be subtracted (because it 
has a negative coefficient) from the total. The number 
of absorbed particles from this term until time t is given 
by —it represents an additional ’’forward 
movement” of the absorbing wall that is ultimately re¬ 
sponsible for the decaying of the random walk in a finite 
time (i.e, avalanches have a maximum value). We end up 
with the problem of an unbiased random walk m{t) with 
an absorbing boundary b{t) moving with the law 

b{t) = -Xmint^^^ + (A8) 

(some non-essential coefficients have been set to 1 in this 
expression). It is straightforward to numerically simulate 
this process, finding the size of the avalanches as the first 
time to at which mito) = &(to); namely when the walk 
is absorbed, and then making statistics on the values of 
to. The results are shown in Fig. [9l Simulations are 
presented for different k values, whereas Xmin are taken 
from its known distribution (Eq. ([SJ). The results give 
a nice verification of the form of the cut off that was 
observed in the full simulations of Fig. [51 

I notice that an analytical solution to the survival prob¬ 
ability of a random walk in the presence of the absorbing 
















FIG. 9: Dotted lines: Size distribution of avalanches described 
by the absorbing of a random walk m{t) by a wall located at 
h{t) = . Results are presented for different 

k values. The values of Xmin are taken from their expected 
distribution in our problem (Eq. ([5])). The continuous line 
is a fitting to the analytical expression (lAlOl) . using a single 
average x^min value. 


wall as in Eq. (jA8l) is not known. However, the form 
of the cutoff can be heuristically determined by inserting 
Eq. (IA8I) in the asymptotic form expected for a solution 
of the diffusion equation of the random walk at large x, 
namely ^ exp(—x^/2t), which gives for the cut off the 
form 


^ exp (-fc^S'^/2 + Xmin/cS') (A9) 

We can thus expect the form of M{S) to be given by 

Af{S) - S-^ exp {-k'^S^/2 + XrmnkS) (AlO) 

where r is given as a function of Xmm in Section II. In 
Fig. [S] we see that this gives a remarkably good fitting 
of the numerical results. From Eq. (lAlOl) it is clear that 
the intensity of the overshoot becomes larger as Xmin 
increases^ Note also that the scaling of Smax as 1/A: is 
immediate form expression (lAlOl) . 
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Eq. (lAlOII seems to indicate that for Xmin = 0 the cut off is the exponent is necessary to fit the results. This indicates 

purely of the form ~ exp(—However numerical sim- that Eq. (lAlOII is certainly not a rigorous formula, 

ulations show that even in this case, a linear in S term in 






